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ABSTRACT 

At fixed stellar mass, the size of low redshift early-type galaxies is found to be a factor of two 
larger than their counterparts at Z'-^l, a result with important implications to galaxy formation 
models. In this paper, we have explored the buildup of the local mass-size relation of elliptical 
galaxies using two visually classified samples. At low redshift we compiled a subsample of 
2,656 elliptical galaxies from the Sloan Digital Sky Survey (SDSS), whereas at higher redshift 
(up to z^l) we extracted a sample of 228 objects from the HST/ACS images of the Great Ob- 
servatories Origins Deep Survey (GOODS). All the galaxies in our study have spectroscopic 
data, allowing us to determine the age and mass of the stellar component. Contrary to previous 
claims in the literature, using the fossil record information contained in the stellar populations 
of our local sample, we do not find any evidence for an age segregation at a given stellar mass 
depending on the size of the galaxies. At a fixed dynamical mass there is only a ^ 9% size 
difference in the two extreme age quartiles of our sample. Consequently, the local evidence 
does not support a scenario whereby the present-day mass-size relation has been progressively 
established via a bottom-up sequence, where older galaxies occupy the lower part this relation, 
remaining in place since their formation. We do not find any age segregation difference in our 
high-z sample either Therefore, we find a trend in size that is insensitive to the age of the 
stellar populations, at least since z~'l. This result supports the idea that the stellar mass-size 
relation is formed at z^l, with all galaxies populating a region which roughly corresponds to 
1/2 of the present size distribution. We have explored two possible scenarios for size growth: 
puffing up or minor merging. The fact that the evolution in size is independent of stellar age, 
together with the absence of an increase in the scatter of the relationship with redshift does 
not support the puffing up mechanism. The observational evidence, however, can not reject at 
this stage the minor merging hypothesis. We have made an estimation of the number of minor 
merger events necessary to bring the high-z galaxies into the local relation compatible with 
the observed size evolution. Since z=0.8, if the merger mass ratio is 1:3 we estimate ~3±1 
minor mergers and if the ratio is 1:10 we obtain ^8 ±2 events. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: formation 
- galaxies: stellar content - galaxies: structure. 



1 INTRODUCTION 

Present-day galaxies show a clear correlation between mass and 
size, with the most massive galaxies having the larger sizes. This 
mass-size relationship has been known both for elliptical and spi- 
ral galaxies for many years. With the advent of large s urveys, like 
the Sloan Digital Sky Survey (SDSS, lYork et al.l2TO0l) it has been 
possible to qua ntify this correlation with high accuracy (see e.g. 
IShen et alj2003 ). However, the mechanisms by which this relation- 
ship is built remain uncertain. For instance, we do not have conclu- 
sive answers to questions like: "were the galaxies bom in-situ at the 
positions where we find them in the local mass-size relation or were 



they bom in another part on this diagram, drifting to their present 
location?". If so, "how much have they grown and what are the 
mechanisms responsible for this displacement?". Answering these 
questions is directly connected to our understanding of how the as- 
sembly of the galaxies has proceeded through cosmic time. In this 
paper we will particularly focus on spheroidal galaxies as it has 
been shown in the last few years that their stellar mass-size relation 
has dramatically changed with redshift. 

Several papers have explored the evolutio n of the stellar mass- 
size relation of spheroid-like galaxies (e.g . iTruiillo et alj|2004l: 
iMcIntosh etai] l2005l : iTmiillo et al] l2006al. |2007| : ISuitrago et al.1 
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I2OO8I : iFerreras et al.ll2009bl ISaracco et al.ll201 lb . In general, they 
all agree with a significant evolution of this relation with red- 
shift. Their results can be summarized as follows: at fixed stel- 
lar mass, spheroid-like g alaxies were significantly more compact 
at higher redshift (e. g. lOaddi et al.' 2005; Trujillo et al. 2006b; 
Longh etti et al] |2007l; Izir metal. 2007; van dg r Wei et al. 2008; 
van D okkum et al. 2008; ICimatti et al.l boOSi ; lOamianov et alj 
20091 ; ,Carrasco et al.„2O10h , with an increase on the effective radii 



by a factor of ~2(4) from z~l(2) (e.g. iTruiillo et al]|2007h . How- 
ever, these observational results say little about the amount of size 
evolution of individual galaxies on the mass-size plane. Neverthe- 
less, at least a few basic statements can be established regarding 
the growth of individual galaxies based on the current observa- 
tional evidence. First, at high-z there are no big spheroidal objects, 
implying that the present-day large elliptical galaxies have either 
formed recently (in-situ) with large sizes or they are the product 
of the evolution of previous compact galaxies that populated the 
high-z stellar mass-size plane. Second, the near ab sence of com- 
pact massive galaxies in the nearb y Universe (Truiil loet alj|200^ ; 
IValentinuzzi et alj2010l ; Pfavlor et al.l2010 ). which were very com- 
mon in the early Universe, indicates that individual objects (at least 
the very old and compact ones) have evolved significantly in size. 

Some recent works have conducted a detailed analysis 
of the buildup of the local spheroid mass-si ze relationship 
jvan der Wei et al.l |2009| ; IValentinuzzi et al.l l2010h . These works 
propose that the formation of this relation is a result of two steps: 
a) the continuous emergence of galaxies as early-type systems 
with larger sizes, as cosmic time increases, due to the de creasing 
availability of gas during their formation phase dtChochfa r & Silk 
l2006l) . and b) their subsequent growth thr ough either gas expul- 
sion in the so-called puffing up scenario jpan et aUl2008l I2OI0I; 
Damjanov et al. 2009) or by minor merging activity ( Naab et all 
2009l ; lHoDkins et alj2009f) . If the above scenario is correct, i.e. that 



the new assembled galaxies are bom with larger sizes as redshift de- 
creases, we should observe that the number density of spheroid-like 
massive galaxies at fixed stellar mass should decrease with increas- 
ing redshift. Furthermore, a gradual change of the age of the galax- 
ies at fixed stellar mass should be expected, in the sense that larger 
galaxies should be younger. However, there is no compelling evi- 
dence of a significant dr op in the number dens ity of elliptical galax- 
ies up to z~l (see e.g. Ferreras et al. 2009b), weakening this for- 
mation scenario. On the o ther hand, in van der Wei et al. ( 2008 ) and 
IValentinuzzi et alj ( I2OI0I) . there is some hint that larger spheroid- 
like galaxies, at fixed dynamical and stellar mass, are younger than 
their compact mass equivalents. 

In this paper we reexamine the buildup of the mass-size rela- 
tionship of spheroidal galaxies with two significant improvements 
in relation to previous work. First, this paper addresses the issue 
of the evolution of early-type galaxies on the size-mass plane by 
comparing a nearby and a distant sample of galaxies, classified and 
analyzed in the same way. We will show in this paper that previous 
studies of the local stellar mass-size relation of early-type galax- 
ies are severely contaminated by galaxies of other morphological 
types. For this reason, the present study is the first one exploring 
objects that have been classified only visually, and not by any other 
criteria, like structural parameters or colours. The second advan- 
tage of the present work is that we have quality spectra for all our 
targets, allowing us - by exploring their spectral energy distribu- 
tions (SEDs) - to obtain reliable star formation histories (SFHs). 
Spectroscopic data is essential to robustly determine the proper- 
ties (stellar mass and age) of the underlying stellar populations in 
both local and high redshift samples, allowing us to make a much 



more consistent assessment of the increase in galaxy size on an in- 
dividual galaxy basis. The information about the ages allows us to 
explore whether the size evolution depends on the properties of the 
stellar populations of these galaxies. This information is key to dis- 
tinguish between the two most likely mechanisms of size growth 
proposed in the literature for elliptical galaxies: the puffing up ver- 
sus the minor merging scenario. We will amply discuss in this paper 
the implications of our findings in relation to these two models. 

The paper is structured as follows. In Section 2 we describe 
the local sample. The connection between stellar age and the mass- 
size relationship is explored in Section 3. In Section 4 we present 
our moderate redshift sample and in Section 5 we quantify the size 
evolution of our galaxies. Section 6 is devoted to explore which 
evolutionary scenario is more plausible according to these obser- 
vations, finally concluding on Section 7 with an overview of our 
results. In this paper we adopt a standard ACDM cosmology, with 
f2™=0.3, Q.K=Q.l and Ho=70 km/s/Mpc. 



2 THE LOCAL SAMPLE 

Our local sample i s taken from the morphological catalogue 
of Nair & AbrahamI (1201 OL hereafter NAIO) obtained via vi- 
sual classification from SDSS imaging. The NAIO catalogue 
com prises 14,034 galaxies fro m SDSS Data Release 4 (DR4, 
Adel man-McCarthv et ai]|2006l) in the redshift range 0.01 < z < 
0.1 with an extinction-corrected g-band magnitude brighter than 
16. From this catalogue, we select the elliptical galaxies (cO, EO 
and E-I-) with a T-Type class -5, resulting in a final sample compris- 
ing 2,656 galaxies with available spectra. 

Our choice for visually classified galaxies aims at minimizing 
the impact of morphological contaminants, which frequently de- 
grade automated classification samples. To illustrate this issue, we 
compare in Fig. [T] the NAIO sample used in t his paper with the 
'Grav es et al.l ( 1200 8l) early-type sample used in Ivan der Wei et al.l 
(200^ for a similar analysis to the one presented here. In con- 
trast with ou r selec tion, which is purely based on morphology, the 
[Graves et al. Illooi) galaxies are selected by the following criteria: 
a) located on the red sequence, b) no emission lines in their spec- 
tra and c) concentration parameter C>2.5. From the comparison 
between both samples shown in Fig. [Tj it follows that a sample of 
early-type galaxies based on the above criteria will be contaminated 
by bulges of both face-on and edge-on spirals. Unfortunately, those 
contaminants are not distributed evenly throughout the sample. In- 
stead, they tend to concentrate in certain regions of the parameter 
space, introducing undesired systematic effects. For example, we 
see that the largest galaxies in the mass-size plane, at a fixed mass, 
are heavily contaminated by spiral galaxies. Another important ad- 
vantage of our selection purely based on morphology is that it is 
neither biased against galaxies with recent star formation activity, 
nor with a passive evolution in their star formation history. This is 
relevant for the purpose of this paper as shown below. We take ad- 
vantage of the fact that virtually all of the galaxies in the sample of 
Graves et al. (2008) have their morpholo gy visually studied by the 
Galaxy Zoo project dUintott et al.ll201 ih . We use their fraction of 
votes for ellipticals (pei); spirals (both clock-wise and anti-clock- 
wise, psp) and edge-on spirals (Pedge) to identify the Late Type 
Galaxy (LTG) contaminants. We have carried out a further visual 
check of these contaminants, finding a complete agreement with 
the results of Galaxy Zoo. Despite the small overall contamination 
rate (~ 1 .8%), the face-on LTGs concentrate in the region with large 



© 2011 RAS, MNRAS 000,[T]-?? 



Dissecting the size evolution of elliptical galaxies since z^l 3 



radii, while tlie edge-on LTGs (--^8%) concentrate towards low val- 
ues of Re and Mdy-n- 

The spectroscopic data and photometric parameters of 
the NAIO sample are retriev ed from the SPSS a rchive. We 
have used spectra from DR7( Ab azaiian et"^ |2009|) . to bene- 
fit from the improved flux cal ibration introduced in DR6 (see 
lAdelman-McCarthv et~ai] |2008|) . The SDSS spectroscopic data 
cover a wavelength range from roughly 3,800 to 9,200 A at an aver- 
age spectral resolution of 3.25 A (FWHM). This instrumental reso- 
lution is not constant but varies in a complex way with wavelength, 
fiber and arrangement. All spectra are both de-redshifted and cor- 
rected for Gal a ctic fo reground extinction, using the dust maps of 
ISchlegel et al] i ll998l) . Hereafter, all size estimates are quoted as 
the circularized effective radius _Re = [h/a)^^^ x RdcV, with 
parameters deVRad_g and deVAB_g taken from the photometric 
SDSS-pipeline. In principle, velocity dispersion data (a) are also 
available from the DRV SDSS-pipeline, although with a moderately 
high ratio of missing values, amounting to over 15% of our SDSS 
sample. Consequently, we have re-calculated the values of velocity 
dispersion with the same spectral fitting method used in this study 
(STARLIGHT, see §3), taking as velocity dispersion the smoothing 
parameter of the stellar po pulation mixture that pro duces the best fit 
to the observed spectrum. iLa Barbera et al. ](|2010^ show that there 
is good agreement between the STARLIGHT and SDSS-DR7 ve- 
locity dispersion values, with only a small systematic trend at the 
low (< 90 km/s) and high (> 280 km/s) ends of the a range. Very 
few measurements (0.4 %) are excluded, with a < 40 km/s, be- 
cause they are considerably smaller than the resolut ion of the base 
SSP models (58 km/s). We have used the J0rgen sen et al.l ( Il995h 
prescriptions to correct the velocity dispersion to the same fraction 
of the effective radius, Re/8, instead of the fixed fiber diameter (3 
arcsec). 



3 THE LOCAL MASS-SIZE PLANE: DISTRIBUTION OF 
GALAXIES ACCORDING TO THEIR STELLAR AGE 



Both lvanderWel et al.l ( l2009h and lValentinuzzi etal.l ( I2OI0I) have 
argued that there is an age gradient within the mass-size plane of 
early-type galaxies: at fixed mass, galaxies with larger sizes are 
found to be younger. As it was argued in § 1 , this observation is ex- 
pected if newly assembled spheroidal galaxies feature larger sizes 
than those systems assembled earlier. For this reason, we have re- 
visited the age distribution of our galaxies in the NAIO sample. 

The age of the stellar populations of our galaxies is esti- 
mated as follows. We use the spectral fitting code STARLIGHT 
jCid Fem andes et aL 2005) to find combinations of single stellar 
population (SSP) models that, broadened with a given velocity dis- 
persion, achieve the best match with the observed galaxy spectrum. 
For the present study, we hav e used the spectral ene rgy distributions 
of the MILES SSP models iVazdekis et Z|l2010l) with a Kroupa 
Universal Initial Mass Function ('Kroupa 200 ly These models 
are based on the MILE^I stellar library ( Sanchez-Blazguez et al.l 
l2006h . which combines both a rather complete coverage of the stel- 
lar atmospheric parameters and a relatively high and nearly con- 
stant spectral resolution, 2.3A (FWHM), optimally suited for the 
spectral resolution of the SDSS data. Our base for the fitting us- 
ing STARLIGHT consists of 138 solar-scaled SSP models with 6 



.iac.es/proyecto/miles 



different metallicities, ranging from Z=l/50 to 1.6 xZ© and 23 dif- 
ferent ages, from 0.08 to 11.22 Gy r. Extinction due to foreground 
dust is modeled with the CCM-law jCardellietal.ll989l) and masks 
are used to avoid emission lines or bad pixels. The Mstars param- 
eter - the fraction of the initial stellar mass which still remains as 
stars at a later time - is extracted from the model predictions and 
used to calculate the stellar mass of our galaxies, Ms. In the present 
work, we have characterized the stellar population mixture of each 
galaxy by its mass- weighted average age, calculated as: 



(age)ii 



where A'^, is the number of SSP models in the base, fij is the 
mass fraction vector, defined as the fractional contribution of the 
SSP with age tj and metallicity Zj, to the total flux, converted into 
mass with the M/Lj of each SSP. Stellar masses, Ms, are also com- 
puted with the fij and M/L_, . Once each SFH has been calculated, 
the corresponding lookback-times are added in order to set all the 
histories to a common z=0 ground. This offset ranges between 0. 1 
(z=0.01) and 1.3 Gyr (z=0.1). 

Figure |2] shows the size-mass correlation of our SDSS sam- 
ple of early-type galaxies, split according to the (mass-weighted) 
ages det ermined by STARL I GHT. We f ollow the spirit of the mod- 
eling of Ivan der Wei et al.l ( I2OO8I) and IValentinuzzi et al.l ( I2OI0I) 
whereby age-segregated samples are expected to occupy different 
regions of the mass-size relation. In order to maximize the differ- 
ence according to age, we only show the upper and lower quar- 
tiles of the age distribution, corresponding, respectively, to galax- 
ies older than 1 1.7 Gyr (black crosses), and younger than 10.2 Gyr 
(grey triangles). These symbols represent the median within bins 
taken at a fixed number of galaxies per bin. The error bar gives 
the RMS scatter within each bin. The figure shows that the pre- 
dicted segregation in the mass-size relation with respect to age is 
not significant when plotted against stellar ma ss (panel a). This 
is in contradiction with the results shown in IValentinuzzi et al.l 
( I2OI0I) (their Fig. 3). We have quantified the separation between 
the young and the old galaxy families by zooming in the re- 
gion where the masses of the galaxies of our two extreme age 
quartiles overlap (Fig. |3}. We quantify the size difference as fol- 
lows: ARe/Re =2(Re,young-Re,old)/(Re,young+Re,old)- As We 

expect from a visual inspection of the figure, the size difference 
between the two families, at a fixed stellar mass, is negligible 
(being compatible with zero change within the statistical uncer- 
tainty). T he reasons for this discr epancy could be several. On the 
one hand. IValentinuzzi et al. 1( 12010) segregates their galaxies using 
luminosity-weighted ages instead of mass-weighted ages as used 
here. Another possibility is that their early-type selection criteria 
based on the automatic code MORPHOT could include a larger 
number of spiral galaxies as contaminants, in contrast with a vi- 
sual classification. We have explored whether using luminosity- 
weighted ages changes our results and find that this is not the 
case. In fact, if we repeat the previous exercise using luminosity- 
weighted ages, we find that the difference between the two extreme 
quartiles is 1.7±0.9% (i.e. very similar to the mass-weighted ages). 
Finally, our ste llar mass-size rela tion is compared with the early- 
type relation o f lShenetal.l ( l2003l) . The agreement is very good for 
objects with stellar mass Ms < 3 x lO^^M©. However, at the 
high mass end we note that t he sizes of our ga laxies are slightly 
larger than those provided b y IShen et ak l (l2003b . a result in agree- 
ment with lGuo et al. I ( I2OO9I) , who find a similar underestimate of 
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Figure 1. Dynamical mass-size relation of the sample used in this paper (NAIO sample; black dots in the right column) compared to the early-type sample 
selection of lGraves et alj flOOSL grey data points). This figure illustrates the different position in the mass-size diagram that spiral galaxies (i.e. contaminants) 
have in this diagram (see text for details). Late-type galaxy contaminants (LTGs) are not distributed homogeneously over the early-type galaxy footprint: 
face-on LTGs mainly live in the top part of the diagram (large radii), whereas edge-on LTGs populate the bottom-left comer (low sizes and masses). 



the Shen et al. sizes of a similar sample of visually inspected early- 
type galaxies. 

Our previous results show that at fixed stellar mass galaxies 
do not show any significant difference in age. However, an inter- 
esting change is found when dynamical masses are considered in- 
stead of stellar ones (panels c and d). In this case, the age seg- 
regation is apparent, with younger galaxies having slightly larger 
sizes. Under the assumption of dynamical homol ogy (i.e. estimat- 
ing t he dynamical masses as Mdyn = ^a^Re/G jCappellari et al.l 
12006) ) the size difference among the two extreme quartiles reaches 
a value of ~16% (~13% in the case of luminosity-weighted ages). 
However, elliptical galaxies are well known for not being an ho- 
mologous family. If we repeat the same analysis using this time the 



dynamical mass accou nting for the non-h omology following the 
expression provided bv lSertin et al.l ( l2002h : 



with 

, 73.32 

K(n) 



+ 0.954, 



(1) 



(2) 



10.465 + (n - 0.95)2 

and n being the S ersic index of the elli ptical galaxies in our sample 
(determined from iBlanton et al.l2005h . we find that the size differ- 
ence, ARe/Re, decreases significantly to (^8% in the case 
of luminosity- weighted ages). Our findings about the size differ- 
ence between the old and young galaxies at a fixe d dynamical mass 
is in q ualitative agreement with the findings by Ivan der Wei et al.l 
l l2009l) . 
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Figure 2. Correlation between effective radius and stellar mass (panel a), 
central velocity dispersion (panel b) and dynamical mass estimated assum- 
ing homology (panel c) and non-homology (panel d) for our local sample. 
The grey triangles (black crosses) represent the youngest (oldest) quartiles 
of the age (inass-weighted) distribution according to our modeling with 
STARLIGHT. The dotted line in panel a) i s the scaling relation of early- 
type galaxies according to lShen et al 
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Figure 3. Zoom into the mass-size distribution of our local sample. The 
regions shown are those where the two extreme age quartiles of the sample 
(young as a dashed line and old as a solid line) have galaxies over a similar 
range of stellar mass (left panel) and dynamical mass (middle and right 
panels). An estimate of the size difference in the overlapping regions are 
given in the top for each panel. 



Given that dynamical mass estimates depend on the velocity 
dispersion quadratically, one would expect that the size difference 
of the galaxies in the two age quartiles could be linked to a change 
on the velocity dispersion between the young and the old popula- 
tion. Panel b) of Fig. [2] confiiTns this point: at fixed velocity dis- 
persion, the older subsample is significantly larger (35.3±0.5%; 
25.3±1.7% in the case of luminosity-weighted ages) than the 
younger galaxies. Alternatively, one could interpret this result as 



follows: at fixed effective radius, older galaxies have lower veloc- 
ity dispersion than their younger counterparts (although the region 
of overlap between old and young galaxies at fixed size is arguably 
rather small). 

Although we agree qualitatively with van de rWeletal.] ( l2009h 
on a size difference between the young and old galaxy families at 
a fixed dynamical mass, our results about the size difference at a 
fixed velocity dispersion is in contrast with them. These authors 
show (in their Fig. 1) that at fixed velocity dispersion the age of 
the galaxies is independent of their size. The reason for this dis- 
crepancy could be double: first, their ages are luminosity-weighted, 
in contrast with our mass-weighted ages, and second, their sample 
suffers from some contamination of spiral galaxies in key places of 
the mass-size diagram. 

Irrespectively of the comparison with other works, our results 
indicate that the size variation due to changes in the stellar pop- 
ulation ages of the elliptical galaxies in the local Universe is very 
small. Although the age trend goes in the direction (i.e. older galax- 
ies being more compact than young ones at fixed stellar mass) that 
one would expect from a progressive bottom-up scenario for the 
buildup of the local mass-size relation model, it is clear that the 
differences in size are very small to be able to reproduce the large 
size variation with cosmic time found at high redshift. We will re- 
turn to this point more extensively in the following sections. We 
conclude that the stellar population ages do not resemble the age 
of the full assembly of the elliptical galaxies and, consequently, 
that after the formation of the bulk of their stellar content, elliptical 
galaxies have experienced a significant evolution in their size. 

3.1 Dynamical structure change of the galaxies with age 

The virial theorem predicts that, at fixed mass, the velocity dis- 
persion will change as the inverse of the root square of the galaxy 
size. Consequently, one would expect, due to the strong size evo- 
lution with redshift observed in the elliptical population, that the 
velocity dispersion of the high redshift objects were significantly 
larger than those found in lo cal galaxies. However, o bservations 
are at odds with this scenario dCenarro & Truiillol2009k : the veloc- 
ity dispersion of the elliptical gal axies, at fixed s tellar m ass, only 
changes moderately with redshift. iHopkins et al. I( l2009bl) have ex- 
plained this mild change in the velocity dispersion suggesting that 
the contribution of the dark matter halo to the gravitational potential 
of the galaxy changes with cosmic time. According to that model, 
in the present Universe the contribution of the dark matter halo on 
settling the velocity dispersion of the galaxies will be higher than 
in the past. 

We can explore whether our local sample shows any hint of a 
dyna mical structure change a s a function of the age as suggested 
by the 'Hopkins et al.l ( l2009b l) idea. To do that we explore both the 
baryonic fraction (top) and the velocity dispersion (bottom) of our 
local galaxies against the age of their stellar populations in Fig. |4] 
The baryonic fraction is defined as the ratio between the stellar 
mass and the dynamical mass, and should roughly correspond to 
the net baryon fraction within the effective radius for our early-type 
galaxies. We bin the sample according to stellar mass as labelled. 
In order to illustrate the effect of non-homology, we show the ho- 
mologous estimates as thick lines and the non-homologous models 
(using dynamical masses from Eq.[TJ as thin lines. 

The strong trend habitually found between the stellar popu- 
lations and velocity dispersion is evident, with the oldest galax- 



ies being the mos t massive ones (see e.g. iBemardi et al 



iGraves et al.l 120081: iRogers et all I2OI0I : iNapolitano et all boioh 



20051: 
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Figure 4. The stellar mass-weighted age is shown with respect to baryon 
fraction (top) or velocity dispersion (bottom) for a range of stellar masses 
as labelled. The points and error bars give the median value and enor in age 
bins chosen at fixed number of galaxies per bin. The thick lines indicate that 
the dynamical masses have been calculated assuming homology, whereas 
the mass estimates assuming non-homology are shown as thin lines. 



wit h a larger dark matter content within the optical radiu s (see 
e.g. lFerreras et aI.l2005l:[Tortora et al.l2009l:lLeier et al.bOllh . This 
trend of an increased dark matter content with galaxy mass is also 
consistent with the results pertaining to whole halos, as shown 
when comparing observed stellar mass functions with cosmolog- 
ical halo abundances (see e.g. Moster et al. 2010). 

We can now probe in more detail Fig.|4] at fixed stellar mass, 
the oldest galaxies feature higher velocity dispersions and a lower 
baryon fraction. The higher velocity dispe rsion for the older galax- 
ies is in agreement with the findings of Cenarro & Trujillo (20091) 
at high redshift and supports the idea of Hopkins et al. (2009b'). 
However, the decrease in the baryonic fraction as a function of age 
seems at odds with the high-z findings. It is interesting to note that 
our estimation of the dynamical mass is done with the sizes found 
in the local Universe, so a direct relation with the baryonic frac- 
tion estimated at high-z is not straightforward. From the analysis 
of the local relation, we find that the age of the most massive lo- 
cal ellipticals is quite homogeneous and also their dynamical struc- 
ture change is limited to Ms/Mdyn'^0.4±0. 1. This suggests that 
the most massive elliptical galaxies formed via an earlier, very ho- 
mogeneous formation process. This scenario is consistent with the 
observed lack of evolut ion in the number density of massive early- 
type galaxies (see e.g. lFontanaetai]|2006l ; iFerreras et ailbOOgbl : 
[Banerji et al. 20^. 

For our intermediate and lower stellar mass bin, ellipti- 
cal galaxies show a much more important trend between age 
and dynamical structure. For instance, we see that for present 
Ms~10^^Mq ellipticals, the baryonic fraction can change between 
0.3 to 0.7 and the velocity dispersion between 150 to 250 km/s. We 
note that our trend, at fixed stellar mass, tow ards a lower baryon 
fraction in older populations is at odds with Ishankar & BemardU 



( |2009|) and iNapolitano eTakl ( |2010|) . They obtain the opposite 
trend, namely m ore dark matter in the younger populations at 
fixed stellar mass dNapolitano et al ]l20T(il) ." or corrected luminosity 
dShankar & BemardilbOOgh T However, our range of stellar masses 
and ages is much shorter, concentrated towards the high mass 
end. Furthermore, the age estimates of Napgiitario_et_al, (2010) are 
based on broadband photometry alone, a method considered robust 
on t he determinations of the stellar M/L but not on age estimates 
(e.g. iFerreras et al ] |2005l) . [Shankar & Bemardi (2009) use instead 
the spectroscopic ages from Gal lazzi et al.l (i2005) who use a combi- 
nation of spectroscopic line strengths. In an independent study car- 
ried out with 40,000 ETGs from SDSS (de la Rosa et al. in prepara- 
tion) sever al methods an d SSP models are compar e d. The method 
of Gallazz i et all ( |2005[) with iBruzual & CharlotI ( |2003|) models 
provide systematically younger ages than the spectral-fitting tech- 
nique with the MILES population synthesis models used for the 
present study. Furthermore, by comparing the performance of the 
model-method combination with repeated observations of the same 
SDSS targets (~2300 repeated spectra), the spectral-fitting ap- 
proach is shown to be considerably more robust than other age dat- 
ing methods. The difference between older and younger galaxies 
may reflect different channels of galaxy formation within the same 
stellar mass bin. As the size of the galaxies is proportional to the 
dynamical mass and inversely proportional to the square of the ve- 
locity dispersion, we obtain, as expected, a slight trend to smaller 
galaxies as a function of the stellar population ages. This trend is 
unable to explain the strong size evolution found at high redshift. 

The detailed analysis of the local mass-size relation reveals 
that the information contained is unable to fully explain which 
mechanisms have followed the elliptical galaxies to reach their 
present sizes. For this reason, it is necessary to conduct a direct 
comparison of the properties of the local galaxies with those of 
equivalent galaxies at high-z to extract such information. This is 
what we do in the following sections. 



4 MODERATE REDSHIFT SAMPLE 

In order to understand in more detail the size evolution of mas- 
sive galaxies and their relation to age, we include in our study a 
sample of visually classified early-type galaxies at moderate red- 
shift (z ;^ 1). The comparison with the local sample allows us 
to probe the evolution of the mass-size relationship over the past 
~8 Gyr The deep image s of the Great Obs ervatories Origins Deep 
Survey (GOODS) fields ( iGiavaUsco et"aLl l2004) taken by the Ad- 
vanced Camera for Surveys (ACS) on board the Hubble Space Tele- 
scope (HST) provide the optimal dataset for visual classification of 
galaxy morphologies out t o redshifts z ^ 1. We use the catalogue 
of early-type galaxies from lFerreras et al1(l2005l) and IFerreras et al.l 
|2009a) in the North and South GOODS fields, comprising 910 vi- 
sually classified early-type galaxies brighter than F775W= 24 mag 
(AB). For a proper comparison with the evolved local sample, we 
need a reliable estimate of stellar age. The broadband photometry 
of the GOODS sample is not good enough for our purposes, and 
we consider a subsample wit h available spectral da ta. The PEARS 
sample of early-type galaxies lFerreras et al.l(l2009ch comprises 228 
galaxies from the GOODS catalogue, with available slitless spec- 
troscopy using grism G800L (HST/ACS). The spectral resolution 
depends on galaxy size, with an average value R= A/ A A ~50 
for our objects. This sample covers a redshift range 0.4<z<1.3. 
The lower redshift was dictated mainly by the requirement of hav- 
ing the 4000A break within the sensitivity range of the grism data. 
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Figure 5. The predicted fraction of stellar mass growth for our local sample 
is shown when backtracked to four redshift bins, as labelled. The histogram 
only considers galaxies with stellar mass > IO'^^Mq. The distribution 
is shown for the youngest (grey, solid line) and oldest (black, dashed line) 
halves, according to the models. 



show the predicted amount of stellar mass for the galaxies in our 
local sample formed since Z'^l according to their SFHs. 

Fig. [5] uses the best fit models from STARLIGHT to quantify 
the net increase in stellar mass from recent phases of star formation. 
We show the mass growth as the ratio between the stellar mass al- 
ready in place at some redshift (Mz) and the current mass at redshift 
zero (Mo) for four redshift bins. The sample is split at the median 
in age measured at zero redshift. One can see that the stellar mass 
growth at z ^ 0.6 stays well below 10% for most of the galaxies, 
especially for the most massive galaxies, which belong mostly to 
the oldest half (black solid lines). 

We have applied to all our galaxies in the local relation the 
evolution in mass predicted from their SFHs and we have rebuilt 
the local stellar mass-size relation taking into account that evolu- 
tion. We consider both the change in stellar age (mass-weighted) 
and the change in stellar mass of the galaxy. For simplicity, we as- 
sume that, within a galaxy, the star formation history does not have 
a radial trend. Our sample does not allow us to probe in detail this 
point, but we note that studies of the colour gradient of early-type 
galaxies at moderate redshift find almost always t he star forma- 
tion concent rated in the centre, i.e. in a blue core dFerreras et al.l 
l2005ll2009ah . The stellar mass-size relation of our PEARS sample 
in comparison with the local sample is shown in Fig. |6] The fig- 
ure shows how the local stellar mass-size relation will look like at 
different redshifts if we correct for the stellar mass evolution. One 
can see that the redshifted local stellar mass-size relation changes 
very little in the high mass regime. The evolution is more evident 
at lower masses, where the galaxies clearly deviate from the local 
relationship. 



In iFerreras et al.l (l2009d) stellar ages are determined using a grid 
of composite models, including chemical enrichment, from which 
best fit ages and metallicities are obtained. However, in order to 
reduce the systematics, we only use this modelling to generate the 
best fit spectra at similar resolution to those from the local sam- 
ple. We note this method should introduce a very small systematic 
given that the values of the reduced obtained for the PEARS 
sample are always of order one, and that the method used in this 
paper to determine ages uses the full SED for fitting, not individual 
absorption lines. 

Ages and stellar masses are re-computed from these spec- 
tra, using the same methodology as f or the local sample (i.e. 
STARLIGHT iCid Femandes et al.il200^ . with the only difference 
being the age range of the model populations. For these galaxies we 
restrict the oldest SSPs to the age of the Universe at the redshift of 
the galaxy. This approach is well justified as STARLIGHT uses the 
full SED to constrain the stellar populations, an equivalent tech- 
nique as the one used with the PEARS dataset. Comparisons be- 
tween STARLIGHT ages and stellar masse s and those determined 
with the chemical enrichment modelling in IFerreras et al.l (l2009cl) 
are fully consistent within error bars. 



4.1 Backtracking the evolution of local early-type galaxies 

By extracting the star formation histories (SFHs) of our local galax- 
ies, one can backtrack their evolutionary paths and estimate the 
amount of new stellar mass created due to the formation of new 
stars as well as the age of the stellar populations at a given redshift. 
To minimize systematic effects, we apply the same methodology to 
both local and distant galaxies to determine their ages. In Fig.|5]we 



5 SIZE EVOLUTION 

We are now in a position to explore the size evolution of the 
early-type galaxies after accounting for the stellar mass growth due 
to new star formation. In fact, the comparison with the observed 
PEARS sample for similar stellar ages will allow us to determine 
the evolution of size at a given stellar mass. 

At the top-left corner of Fig.|6] the local sample is shown using 
the same criterion as in Fig. |2] with individual galaxies shown as 
small dots. We include in that p anel the loc al trend of SDSS early- 
type galaxies (long dashed line. iShen et al.l i,2003), showing agree- 
ment with our local sample, except for the most massive end, as 
discussed in Section 3. In the following panels, PEARS individual 
galaxies appear as solid (open) circles, with ages younger (older) 
than the median within each redshift bin. The standard downsizing 
trend is apparent in this figure, with the younger PEARS galax- 
ies having the lowest s tellar masses. If the proposed model in 
Ivan der Wei et al.l ( |2008|) were correct, with the youngest galaxies 
being more extended, at a given stellar mass, than the older coun- 
terparts, one would expect this segregation to be more evident at 
higher redshifts, where the effect of lookback time makes it easier 
to discriminate with respect to age (i.e. a reduced age-metallicity 
degeneracy). However, no clear trend with respect to galaxy size is 
found in our data. 

Our best fit models for the local sample predict very small 
stellar mass changes (see Fig. |5]l, at levels that correspond to 
A logMs 0.05 dex along the horizontal direction in Fig. |6| The 
comparison with the PEARS sample shows that there is a notice- 
able "vertical" evolution (i.e. change in size). This one can be il- 
lustrated by comparing the (redshift zero) size of the local galaxies 
with the observed size of the PEARS galaxies, within subsamples 
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Figure 6. Comparison of the stellar mass-size relation for the extrapolated values of the local sample and the PEARS sample, in redshift bins, as labelled. The 
black crosses (grey tiiangles) represent the oldest (youngest) half in mass-weighted age of the local sample extrapolated backwards in time according to the 
best fit star formation history (see text for details). The circles correspond to the PEARS sample. Open (sohd) circles represent galaxies older (younger) than 
the median value, computed within each redshift bin. Note the local sample is corrected regarding the evolution of the total stellar mass, but not in size, so that 
a vertical shift should be expected when compar ing both samples (see figure [T]. At z=0 (top-left panel) individual galaxies are shown as dots and the local 
early-type stellar mass-size relation from SDSS dShen et al.l2003l) is given as a long-dashed line in all the panels. 



of the same stellar age. Fig. [7] shows the size evolution for galax- 
ies with stellar mass in the range 5 x lO"' <Ma/M0 < 3 x lO". 
We have fitted our evolution using the following parametrization: 
R(z)=R(0)(l + 72) wit h R(0) the size ob tained from the local stel- 
lar mass-size relation jShen et al.ll2003l) . Our data is compatible 
with the following value 7=-0.657±0.122 for the full galaxy sam- 
ple and with 7=-0.631±0.176 for the young subsample and 7=- 
0.674±0.160 for the older subsample (uncertainties quoted at the 
68% confidence level). 

Fig. |7] sho ws that the si z e evo lution is significant, in agree- 
ment with, e.g. iTruiillo et all ( l2007t) . with galaxies at 2 ~ 1 be- 
ing ~50% smaller in size than their local counterparts. Notice the 
little difference between the trend of the sample segregated with 
respect to age (large open/solid grey circles). This is one of the 
most important result of this work and implies that the amount of 
size evolution that elliptical galaxies suffers since z~ 1 is indepen- 
dent of the age of the galaxies at each redshift interval. This means 
that, the full population of elliptical galaxies, independently of its 
level of star formation, experiences a similar evolutionary mecha- 
nism for assembly. This is once more a result in contradiction with 



the idea that younger galaxies at all redshifts are bom with signif- 
icantly larger sizes than their older massive counterparts . In other 
words, our results point out to a similar displacement in the stellar 
mass-size relation of all the galaxies in the sample (independently 
of their age). 



6 CONSTRAINING THE DIFFERENT EVOLUTIONARY 
PATHS OF THE ELLIPTICAL GALAXIES SINCE Z-1 

In this section we explore the current most likely scenarios pro- 
posed to explain the evolution of elliptical galaxies on the mass- 
size plane. We use the results obtained here and in previous papers 
to constrain those scenarios. 

In what follows we consider that both the size and the stellar 
mass growth of the elliptical galaxies can be described as the con- 
tribution of three different processes: i) formation of new stars in 
the galaxies as a result of gas consumption, ii) accretion of already 
formed stars from merging of different subunits and iii) gas ejec- 
tion from the activity of either an AGN and/or supernova galactic 
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Figure 7. Size evolution of the PEARS galaxies split according to their stel- 
lar age. Only galaxies with stellar mass between 5 X 10^" and 3 X IO^Mq 
are considered. Small open (solid) circles show the evolution for galaxies 
older (younger) than the median within each redshift bin. The diamonds 
give the size evolut ion of massive (Mj, > IO^^Mq) spheroids (Sersic in- 
dex n> 2.5) from iTruiillo et al\ bOOTi) . The error bars of the value of 
iTruiillo et al.l j2007l) represent the scatter of their sample. The lines rep- 
resent a hnear fit R(z)=R(0)(l + 72) to the different galaxy populations: 
solid line (older subsample) and dashed line (younger subsample).The grey 
solid area is the fit to the full population including the 68% confidence level. 



winds. We parameterize the effect of these three processes in the 
mass and size of the galaxies as follows; 

AMs = AMs,SF + AMs.acc (3) 
Are = Are.SF + Are.acc + Are,agn (4) 

with AMs,sF and AAIs.acc representing the increase of the stellar 
mass due to star formation and by accretion of new stars into the 
galaxies, respectively. Are,sF, Are,acc and Ar^^agn correspond to 
the increase in size by star formation, accreted stars and by expan- 
sion due to galactic winds either created by the effect of a central 
AGN or supemovae explosions. 



6.2 Puffing up model: AGN and/or supernova galactic winds 
effects 

iFan et al] ( l2008l |2010!) have proposed a mechanism based on the 
removal of gas as result of AGN activity to explain the size growth 
of early-type galaxies. According to these authors, the rapid ex- 
pulsion of large amounts of gas by quasar winds destabilizes the 
galaxy structure in the inner, baryon-dominated regions, and leads 
to a more expanded stellar distribution. A similar idea - but based 
on the ga s expulsion associated t o stellar evolution - has been pro- 
posed bv lPamianov et af] ^20091) . The prediction from the puffing 
up model can be parameterized as follows: 



AMs,sF = (5) 

AMs,acc = (6) 

AMs = (7) 

Are.SF = (8) 

Ar,,acc = (9) 

Are = Are.agn (10) 



In other words, all the galaxies in the stellar mass-size relation 
should just evolve vertically in this relation without any increase of 
stellar mass. Consequently, the size evolution we observe at fixed 
stellar mass should be directly interpreted as the total size evolution 
of the galaxies. 

This model agrees with observations at predicting a small for- 
mation of new stars due to the removal of gas from the galaxies. In 
addition, this model fits well with the lack of evidence of significant 
evolution in the number density of massive ellipticals since z~I. 
However, we find that our data is in conflict w i th the model in sev- 
eral aspects. First, according to the Ipan et al.l (l2008h model, after 
the formation of the compact structure, the AGN activity will re- 
move the gas, triggering a fast growth process (~20-30 Myr based 
on recent simulation^: iRagone-Figueroa & Granatd ( l20Ilb . This 
would imply that galaxies with stellar populations older than ~1 
Gyr should be already located in the local stellar mass-size rela- 
tion. This is not what our data shows. We have galaxies (old and 
young) at the same distance from the local relation at all redshifts. 
For instance, at z=I the mass-weighted age of our sample is 3.9 Gyr 
for the old subsample and 3.5 Gyr for the young subsample. We can 
consequently assure that the mechanism that is operating in the size 
evolution of our galaxies does not know about the age of the stel- 
lar populations. This is in contradiction with the puffing up model. 
In addition, a natural prediction from the puffing up model is that 
the scatter of the stel lar mass-size relation will increase with red- 
shift (Fan et alj|2o"lol) . with some galaxies already in place on the 
local relation and others still in a very compact phase. We do not 
observe any increase in the scatter of the stellar mass-size relation 
with redshift in our data. 



6.1 Observational facts 

The results of this paper show that AMs,sf is very small (i.e. 
AMs^SF « Ms) and also that the evolution of the size of the 
galaxies is quite independent of the age of their stellar population, 
so Are (old) ~ Are (young). Due to the little increase in the stel- 
lar mass due to in-situ star formation we assume from now on, to 
simplify the discussion, that, if any, AMs ~ AMs.acc for the el- 
liptical galaxies since z~ 1 . 



6.3 Major dry mergers 

Major mergers (i.e. mergers of galaxies with similar mass) 
were first considered as one the likely paths for size growth 
in elliptical galaxies. Major dry mergers can increase the 
size in a wa y almost directly propor ti onal to the mass 
increase (e.g. Iciotti & van Albada I l2()0ll: iNipoti et al.l I2OI0I: 

^ In the case of supernova winds an important mass loss event could last 
even ~0.5-l Gyr. 
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iBovlan-Kolchin et alj20o3 : lNaab et alj2007l') . This evolution is not 
strong enough to be compatible wi th the low number of major 
mergers observed at least s i nce Z'^X jBundv et al.l200^ : IWild et aP 
1200^ ; Ide Ravel et alJIlOOsl i lLopez-Saniuan et al.ll20ld). a s well as 
with recent numerical simulations ("Khochfar & Silkl2009t) . For this 
reason, we will not consider this mechanism further. 



6.4 Minor dry merging 

Another possible scenario for elliptical gal axy growth involves 
to m i nor mergers on p a rabolic orbits (e .g. Khochfar & Burkerg 
|2006| ; iMaller et aLlbood iNaab et alj|2009l : iHopkins et alj|2009bh 
Through this channel, the new accreted stars as well as the redistri- 
bution of stars in the main galaxy, preferentially populate the outer 
region of the objects. For this reason, this me chanism h a s been con- 
sidered a_verY_efficieiUway of size growth. lFanetal.l ( l2oTol) . fol- 
lowing |Na^^ ai]( l2009l) . show that the fractional variation of the 
gravitational radius and the velocity dispersion of the main galaxy 
before (i) and after (f) a minor merger is: 



R^ 



1 +772-' 
1+77 



(11) 

(12) 



with T] defined as M/=Mi(l+?7) and a representing the exponent 
of the local stellar mass-size relation (R=bM"). Shen et al, . (2002) 
proposes a faO.56 with b=2.88x 10""+"° (in units of 10"Mq). 
On what follows, we implicitly assume that the gravitational radius 
is proportional to the effective radius of our galaxies. This is only 
strictly correct as long as the galaxies do not change the shape of 
their surface brightness profiles during the minor merger process. 

It can be shown that after N mergers of equal mass ratio 77, the 
final mass, velocity dispersion and radius can be written as: 



Rl 



M 



1 + 77'- 



i + n 



(13) 
(14) 

(15) 



We can now make an estimation of the number, N, of minor merg- 
ers a galaxy requires in order to reach the present stellar mass-size 
relation. The final size of the galaxy can be written in terms of the 
initial size, the size evolution at a fixed stellar mass (provided by 
the observations) and the difference in stellar mass as follows: 



log Rf = log i?i + A log R 



+ a\og{Mf/M^). 



(16) 



The evolution at fixed stellar mass for different redshifts is deter- 
mined by the size evolution found in our data A log R\m^ f^^^^ ~ 
~ log(l + 72), with 7 = -0.657 ± 0.122. Using Eq. (13)7(14) 
and (16) we find for the number of minor mergers: 



N = 



-log(l + 72) 



log 



1+1?^-° 



(17) 



We show in Fig. [8] this number as a function of redshift for two 
different values of r\: 1/3 and 1/10. As expected, the number of mi- 
nor mergers is a function both of redshift and the mass increase per 
merger, t). We can use these estimations in the number of minor 




Figure 8. Predicted number of minor mergers as a function of redshift ac- 
cording to the observed size evolution in our data. The grey regions repre- 
sent the number of mergers within a 1 ct uncertainty, for two choices of the 
mass ratio (77). As expected, the number of minor mergers is a function of 
r\, being larger for smaller ratios. 



mergers as a function of redshift to determine the increase in size, 
stellar mass and velocity dispersion that individual galaxies suffer if 
their evolution is dictated by the minor dry merging scenario. This 
is quantified in Table [T] Since z~0.8, individual objects undergo 
size growth by a factor around 3.5, whereas the stellar mass grows 
a factor around 2.5. As expected, the velocity dispersion of the in- 
dividual galaxies decreases with time due to the minor merging. 
The evolution is, however, very mild. We can compare this evolu - 
tion with the observed values found in ICenarro & Tniiillol (l2009l) . 
The comparison is not straightforward as those authors measure 
the velocity dispersion increase at a fixed stellar mass, whereas we 
have followed the evolution for individual gala xies and we observe 
that th e increase in stellar mass is not negligible. ICenarro & TruiiUol 
( l2009h find a velocity dispersion change ranging from 0.84 to 0.90 
since z=0.8 which is in good correspondence to the observed values 
predicted here from the observed size evolution. 



7 DISCUSSION AND CONCLUSIONS 

In this paper we have explored how the local stellar mass-size re- 
lation of elliptical galaxies has been built up since z~I. We have 
compiled a sample of visually classified elliptical galaxies since 
z~l from the GOODS datasets, as well as from SDSS data. All 
our galaxies have spectroscopic data that enable a robust constraint 
of the age and mass of the underlying stellar populations. Both the 
study of the fossil record in the local relation as well as the analysis 
of the stellar mass-size relation evolution with redshift agree on an 
evolutionary mechanism that is mostly insensitive to the age of the 
stellar populations of the galaxies at all redshifts. 

We do not find any clear evidence for a progressive buildup of 
the local stellar mass-size relation following a bottom-up sequence. 
In other words, we do not observe that the smaller galaxies, at fixed 
stellar mass, are generally older than the larger galaxies. On the 
contrary, the local stellar mass-size relation seems to be in place 
(with a similar slope and scatter, at least since z~l) but with all the 
galaxies presenting a "vertical drift" towards smaller sizes. 

The analysis of our data rejects the puffing up scenario which 
proposes that the growth in size is due to the rapid expulsion of 
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Table 1. Number of minor mergers, size, mass and velocity dispersion evolution for individual galaxies according to the minor dry merging model compatible 
with our data. 



Redshift 


<N>(±l(j) 


R//R,(±1(t) 


M/ /Mi (±1(7) 


(T//(Ti(±l(T) 


77=1/3 


0.2 


0.62±0.12 


1.27±0.06 


1.20±0.04 


0.97±0.01 


0.4 


1.35±0.29 


1.70±0.19 


1.48±0.12 


0.93±0.01 


0.6 


2.24±0.53 


2.44±0.50 


1.93±0.29 


0.89±0.02 


0.8 


3.38±0.92 


3.95±1.35 


2.73±0.71 


0.84±0.04 


1.0 


5.00±1.64 


8.44±4.74 


4.70±2.06 


0.78±0.06 


»7=1/10 


0.2 


1.40±0.28 


1.24±0.05 


1.14±0.03 


0.96±0.01 


0.4 


3.02±0.65 


1.61±0.16 


1.33±0.08 


0.91±0.02 


0.6 


5.00±1.20 


2.21 ±0.40 


1.62±0.18 


0.86±0.03 


0.8 


7.56±2.05 


3.39±1.04 


2.09±0.40 


0.80±0.05 


1.0 


11.20±3.66 


6.61±3.39 


3.09±1.04 


0.72±0.08 



large amounts of gas by the effect of an AGN or supemovae-driven 
winds. In fact, two key predictions of this model, the increase of 
the scatter in size with redshift at fixed mass, and the absence of 
old galaxies with small sizes are not observed in our data. Our 
data, however, is not in conflict with an increase of the galaxy sizes 
through minor merging. Minor merging has been also favoured by 
studies us ing different methods and with samples at hig her red- 
shift (i.e. iBezanson et al.|[2009l : Ivan Dokkum et alj 12010) . Under 
this hypothesis we have calculated the number of minor mergers 
that would be necessary to build the local stellar mass-size relation 
in agreement with the observed size evolution. Since z=0.8, we find 
~3±1 mergers with ratio 1:3 or ~8±2 with ratio 1:10. 

The data analyzed in this work together with the evidence 
collected in recent papers (e.g. iKavirai et al.ll200^ : IShankar et al.l 
[2OIO; Nierenberg et al. 201^ only leaves the minor merging sce- 
nario as a viable mechanism for the size increase of elliptical galax- 
ies at least since z~l. Proving ultimately, however, that elliptical 
galaxies grow by minor merging will require a direct quantification 
on the minor merger events found in high redshift galaxies and an 
exploration of the age and metallicity gradients of the stellar popu- 
lation in local elliptical galaxies. 
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